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Fermi level pinning at the Ge(OOl) surface — A case for non-standard 
explanation 
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To explore the origin of the Fermi level pinning in germanium we investigate the Ge(OOl) and Ge(001):H 
surfaces. The absence of relevant surface states in the case of Ge(001):H should unpin the surface Fermi level. 
This is not observed. For samples with donors as majority dopants the surface Fermi level appears close 
to the top of the valence band regardless of the surface structure. Surprisingly, for the passivated surface 
it is located below the top of the valence band allowing scanning tunneling microscopy imaging within the 
band gap. We argue that the well known electronic mechanism behind band bending does not apply and a 
more complicated scenario involving ionic degrees of freedom is therefore necessary. Experimental techniques 
involve four point probe electric current measurements, scanning tunneling microscopy and spectroscopy. 
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I. INTRODUCTION 

The renewed interest in germanium is driven mostly 
by the prospects it offers to microelectronicsi*^. To this 
end a detailed understanding of surface electronic and 
structural properties is essential. The position of the 
surface Fermi level (SFL) and the structure of surface 
states are key features in this context. 

There is a generally accepted theory of the Fermi level 
pinning (FTP) in semiconductors, formulated decades 
ago^, which operates within the electronic degrees of 
freedom only. This understanding is useful in explain¬ 
ing many phenomena, e.g independence of the SFL from 
the bulk doping and creation of two-dimensional electron 
gases^. For germanium surfaces this model was adapted 
by P. Tsipas et al. in Ref. Nevertheless, the strength 
and nature of FLP in germanium is still under debate as 
outlined recently^. 

In our recent paper— we demonstrated that the 
surface Fermi level is pinned close to the valence band 
maximum for Ge(OOl). In that letter we concentrated 
on the conductance data and their interpretation. We 
suggested that the SFL is pinned by the surface states 
(dangling bonds) due to the fact that there are no addi¬ 
tional structures on the surface capable of pinning the 
SFL. Our conductivity experiments reported below are 
consistent with the surface Fermi position revealed by 
Scanning Tunneling Spectroscopy (STS). Surprisingly, 
these data suggest a metallic character for the Ge(001);H 
surface resulting from the SFL pinned below the top of 
the valence band. What we find of particular interest is 
that the surface states recovered from STS seem to be 


^lElectronic mail: j.lis@uj.edu.pl 


at odds with the well-known mechanism behind FLP. 
We also review experimental reports, mostly based on 
angle resolved photoelectron spectroscopy (ARPES), 
suggesting a failure of the electronic theory of the FLP 
for germanium. We postulate that a mechanism similar 
to that valid for silicon is in effect. 


II. EXPERIMENTAL 

The experiments were performed on three different 
germanium (MTI Gorporation) samples: undoped 
(~ 45 Ul cm), n-type Sb-doped (0.008 — 0.009 Ul cm) and 
p-type Ga-doped (0.1 — 0.5 Ul cm). The intrinsic doping 
for germanium is n-type, with the density of carriers 
about 10^^ cm~^. The experiments were done in an 
ultra-high vacuum (UHV) system with a base pressure 
below 5 ■ 10“^° mbar. The sample preparation consisted 
of subsequent sputtering (Ar+, 800 eV) and annealing 
to 1040 K cycles. The cleanness of the (001) surface 
was checked with the low energy electron diffraction and 
the scanning tunneling microscopy (STM). Afterwards, 
to passivate the surface, the substrate was exposed to 
atomic hydrogen flux dosed with a home-built hydrogen 
cracker (10“^ mbar partial pressure, sample temperature 
480 K). The procedure resulted in a well-hydrogenated 
surface with the defect density down to 2-4% of surface 
dimers, which translates to a defect density up to 4-10^^ 
cm“^. Further details of the preparation process may 
be found in Ref. H. Prepared samples were trans¬ 
ported in the UHV conditions to the surface transport 
measurement chamber or to the STM chamber. The 
surface transport measurements were carried out at 
room temperature in UHV Nanoprobe System (Omicron 
Nanotechnology GmbH) equipped with a high resolution 
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scanning electron microscope (Gemini Column, Carl 
Zeiss). Scanning tunneling spectroscopy and microscopy 
experiments were done at the liquid helium temperature 
(4.5 K). 


III. EXPERIMENTS AND DISCUSSION 
A. Resistance analysis 

In the previous paper— we reported resistance measure¬ 
ments on clean Ge(OOl) for three samples: n-type, p-type 
and intrinsically doped. The results were divided into 
two qualitatively different groups: for n-type samples two 
dimensional current flow was observed while for p-type 
samples a three dimensional current character was seen. 
For the explored distances between the electrodes the 
current confinement in a thin layer beneath the surface 
is not a trivial effect resulting from finite sample width 
as its occurrence correlates with the doping type. The 
reduced dimensionality of the current was interpreted as 
a result of FTP followed by the creation of an inversion 
layer beneath the surfaces on n-type doped samples. 

In that paper— we discussed the dimensional assignment 
in more details. Starting from the current flow equation, 
we showed that current confinement in the direction nor¬ 
mal to the surface results in a scale independence of re¬ 
sistance measurements. In the case of three-dimensional 
currents the results depend on the absolute distances be¬ 
tween electrodes. The experimental points closely repro¬ 
duced the relevant scaling laws^. 

Here we report on resistance measurements for 
Ge(001):H in the same experimental settings as in Ref. 0. 
We applied the four point probe method. The current 
source and drain were placed at distance D (between 2 
and 16 //m) while the distance between two probes mea¬ 
suring the voltage drop was changed. All measurements 
were done in a colinear geometry symmetric with the 
midpoint of the segment set by the current supplement¬ 
ing electrodes. In these settings the formula for mea¬ 
sured resistance in the case of two-dimensional current 
flow yields 



FIG. 1. Data acquired on n-type doped sample for pure 
Ge(OOl) and after surface passivation for the distance D = 
8 fim between the current source and drain. The distance x 
on the horizontal axis shows the distance between electrodes 
measuring the voltage drop as a fraction of D, see Ref. 0. For 
different surfaces, the measurement outcomes do not differ 
much. Fitting to formula o gives two-dimensional conduc¬ 
tivity cr 2 = (2.61 ± 0.03) • 10“^n/r2 for the bare surface and 
CT 2 = (2.57 ± 0.02) • 10“'^n/fl, which is close to the values 
reported in Ref. 0 The smooth curve shows the function o 
with the parameter CT 2 fitted to the Ge(001):H data. 

a 3 are expressed using different units. 

Fig. [T] and [5] show representative data obtained for the 
Ge(OOl) and Ge(001):H surfaces for n-type and p-type 
doped samples (similar agreement is observed for the in¬ 
trinsic sample). There are small variations of the conduc¬ 
tance introduced by surface treatment. These differences 
are within conductivity fluctuations observed for samples 
cut from one wafer. Consequently, the dimensional char¬ 
acter of the current flow is not altered. This provides a 
hint that the SFL is insensitive to the actual structure of 
surface states. In other words, it suggests that the inver¬ 
sion layer stays intact upon hydrogen termination of the 
surface. 


B. Scanning tunneling spectroscopy and microscopy data 
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where <T 2 is the two-dimensional conductance and x is the 
distance between the probing electrodes normalized with 
D. For the three-dimensional currents the formula reads 


R = 
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where as stands for the three-dimensional conductivity. 
A more detailed analysis of the above relation may be 
found in Ref. 0. These formulas highlight the qualita¬ 
tive difference between the two-dimensional and three- 
dimensional current flows, suffice it to say that a 2 and 


We concentrate on data acquired for the intrinsically 
doped sample. As it is effectively n-type doped, the bulk 
Fermi level is located above the mid-gap. For n-type 
doped sample the results are very similar. The case of 
the p-type doped sample is not that interesting, as the 
band bending is less obvious. The STS data obtained 
for bare Ge(OOI) and hydrogen passivated surfaces are 
shown in Fig. [31 In both cases there is a plateau corre¬ 
sponding to the surface band gap. It is about 0.2 V wide 
for Ge(OOI) and 0.85 V wide for Ge(00I):H, see inset in 
Fig. |0and Ref. 0 . For the passivated surface dl/dV for 
the zero bias voltage does not vanish. It suggests that the 
SFL is located below the top of the valence band. Indeed, 
we can repetitively acquire STM images for any positive 
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FIG. 2. Representative data acquired for p-type doped sam¬ 
ple for pure Ge(OOl) and passivated surface for D = 8 /rm. 
For different surfaces, the measurement outcomes are sim¬ 
ilar, lower resistance for Ge(OOl) than for Ge(001):H is 
not systematically observed. Here, fitting results in ua = 
(1.62 ± for the bare surface and 0-3 = (1.38 ± 

for the hydrogenated one. The smooth curve 
shows relation (l2| with the parameter fitted for Ge(001):H. 
Poor htting quality for a; ~ 1 can be explained^ as a result of 
conductivity variation near the surface. 



FIG. 3. Spectroscopic data for Ge(001)-p(2x2), Ge(OOl)- 
c(4x2) and Ge(001):H acquired at 4.5 K, the feedback set- 
point in all cases was -0.5 V and 500 pA. The arrows mark 
the onset of the valence band. Note that the signal is ex¬ 
tremely weak, unnoticeable when bands are observed. The 
data in a broader range showing the surface band gap for the 
passivated surface can be found in Ref. [^. The quick change 
of dl/dV visible in the inset for the bare surfaces for about 
0.4 V starts about 0.2 V in the close view. 


bias voltage at the Ge(001):H surface, see Fig. H) The 
resulting pattern coincides with the contrast obtained at 
negative biases. This makes us believe that the SFL in 
this case is located slightly (10-30 meV) below the top of 
the valence band. 

In the case of Ge(OOl), the SFL is also located at the top 
of the valence band, as suggested in Fig. 3. We cannot 
obtain repetitively STM images for small positive biases 
(up to 0.2 V). However, on rare occasions we happen to 
arrive at tip apexes that admit such STM imaging. This 
suggests that STS and STM experiments allow determi¬ 
nation of the SFL within several tens of milivolts accu¬ 
racy. Such indeterminacy appears to be an intrinsic fea¬ 
ture of the method linked to tip-induced band bendingi^. 
Furthermore, the bare surface is much more chemically 
reactive than its passivated counterpart. As such, the 
tip probably cannot approach the bare surface as close 
as the passivated surface, which may also contribute to 
the difference between the SFL observed in both cases. 

To summarize, the SFL seems to be pinned at the top 
of the valence band nearly unaltered by the passivation 
of Ge(OOl). This is in agreement with photoelectron ob- 
servationsii. We remark that the insensitivity of the SFL 
to passivation explains conductance measurements in the 
same way as in Ref. Indeed, the inversion layer is not 
destroyed by hydrogen adsorption. 

C. Surface state structure 

There are two low-temperature Ge(OOl) reconstruc¬ 
tions: p(2x2) and c(4x2). In both cases the surface 



FIG. 4. STM images of Ge(001):H acquired for -0.2 V bias 
voltage and - 1 - 0.2 bias voltage, upper and lower panel respec¬ 
tively. The nearly indiscernible contrast appears as a result 
of imaging the valence band states above the SFL. 

atoms reduce the number of unsaturated bonds by form¬ 
ing bonded pairsi^. These dimers form long rows. Two 
atoms in a pair are not equivalent which is manifest in 
a tilted geometry. The rows consist of alternately tilted 
pairs. The two reconstructions result from two possible 
relative orientations of pairs in adjacent rows. These sur¬ 
face arrangements appear nearly equivalent energetically 
and can be reversibly switched in STM experiments^^. 
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Both experimental and theoretical findings regarding the 
Ge(OOl) surface point to the following structure of surface 
bands^. The lower surface band, called usually D^p, is 
located very close to the top of the bulk valence band and 
it is a surface resonance. The name D^p points to the fact 
that the related electron density is concentrated on the 
upper dimer atoms. This band splits into two branches 
when going from T to J 2 for either reconstruction. The 
electric neutrality of the surface requires that this band 
be fully occupied. There is an upper surface band about 
0.3-0.45 eV above the lower band. It is termed Ddown 
due to maximal electronic density concentrated on the 
lower atoms in the dimers. This band stays in the bulk 
band-gap, which for germanium reads 0.742 eV at 0 K 
and 0.661 eV at room temperature. We note that in var¬ 
ious studies the SFL always appears between these two 
surface bands at variable positions. 

The data shown in Fig. 3 correlate well with the above 
description. The surface vs. bulk character of the occu¬ 
pied states just below the SFL cannot be judged from our 
experiments. However, about 0.4 eV above SFL there is 
a peak that we interpret as the surface band since it is 
located within the bulk band gap at a position where 
a surface band is found in the ARPES experiment s^ ^ 1 
and is washed away upon passivation. 

The only feature appearing in STS data we cannot firmly 
identify is a very weak and flat peak observed just above 
the SFL for Ge(001)-c(4 x 2). It is worth noting that 
one can find such a structure in the STS data published 
by Gurlui^. It seems unlikely to be a surface band. It 
appears right above the valence band, where shallow ac¬ 
ceptor states are located. The shallow elemental accep¬ 
tors have an energy about 10-15 meV above the edge of 
the valence band, two vacancy ionization energies are lo¬ 
cated 0.05 eV and 0.14 eV above the valence bandit and 
for the vacancy-Sb complex transition states are located 
between 0.11-0.18 above the top of the valence band^. 
These energies correlate with the peak width. However, 
it is not clear why those states are evident for this par¬ 
ticular surface reconstruction only. 

Hydrogen termination of the Ge(OOl) surface does not 
alter bonds between the surface Ge atoms and, conse¬ 
quently, does not touch the rows on the surface. One 
hydrogen atom binds to every surface Ge atom turning 
the dimer atoms equivalent. The hydrogen terminated 
Ge(OOl) surface has not attracted as much attention as 
the clear surface, but there are several studies available 
such as both ARPES^i results and recently published 
STS dataii^. They point out similar surface state struc¬ 
ture as in Fig. 3. The surface band gap increases to 
about 0.85 eV. ARPES data show clearly that the occu¬ 
pied surface states are located deeply below the top of 
the valence band (about 3-4 eV). There is no additional 
surface-related feature close to the Fermi level seen in the 
ARPES experiments. The unoccupied surface bands are 
above the lower bulk conductance bands^. 

The above resume shows that the spectroscopic exper¬ 
iments explore the actual structure of the surface states. 


The absence of any surface state for Ge(001):H within the 
bulk band-gap makes the observed band bending a puz¬ 
zle. Noteworthy, the results obtained in the four point 
probe experiments at Ge(OOl) and Ge(001):H show that 
dangling bonds (or any surface channels) do not consid¬ 
erably contribute to the electron transport even at sub¬ 
micron scale. 


D. Fermi level pinning 

As already mentioned, band bending at surfaces is due 
to the presence of the surface states not aligned with the 
bulk bands. In general, the bulk position of the Fermi 
level leads to high occupation of the surface states and, 
consequently, the surface gets charged. As a result, an 
electrostatic field appears changing the position of the 
surface Fermi level. In turn, the positions of bands are 
shifted with respect to the Fermi level at the surface and 
the surface charge is significantly reduced. 

Having said this, it is evident that there is no reason 
for the Fermi level to appear at the top of the valence 
band on Ge(OOl). For the p-type doped sample, this 
could happen for high enough doping. But for the n-type 
samples, the bulk Fermi level must be above the mid-gap. 
At T=5 K the Fermi-Dirac factor excludes any significant 
thermal occupation of states located E=13 meV away 
from it, as 


exp 


E 



7.4 • lO"^'^ 


while the surface band density is of the order of 6 • 10^^ 
states per cm^. Hence, any occupation of the empty sur¬ 
face band on Ge(OOl) can be excluded. Furthermore, 
for SFL positioned 0.2 eV higher this occupation would 
not be significantly different. In this context, we ana¬ 
lyze three possibilities. First, for the neutral surface the 
mechanism behind sticking of the SFL to the top of the 
valence band is not clear. Indeed, the theory of FLP 
gives rise to the following formula (Schottky approxima¬ 
tion^) relating the surface charge density Sg needed for 
band bending (j) and bulk dopant density Ni, 


5s ^ \/ (jiNb- 


(3) 


The formula is valid within the band gap. For the neu¬ 
tral surface i5s = 0 the band bending also vanishes. But 
experimentally it yields at least (f) « 0.3 —0.4 eV. Second, 
one can think of superficial surface states located below 
the top of the valence band that, when occupied, make 
the surface negatively charged. But there are good rea¬ 
sons to reject the presence of such states as neither theory 
nor experiment hint at such a possibility. These states 
should be present both on the bare and passivated sur¬ 
face. However no surface states are observed in ARPES^i 
of the hydrogen terminated surface. Third, the surface 
can be positively charged due to the states from the va¬ 
lence band above the SEL. But positive charges do bend 
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bands downwards to move the SFL away from the valence 
band. In the case of n-type doped germanium upward 
shifting of the bands is evident. 

Using eq. m we estimate that for the intrinsic slab the 
charge density needed to cause band bending of the order 
of 0.3-0.4 eV yields 5 • 10® cm“®, and for the doped sam¬ 
ple this is 4 • 10^^ cm“®. These estimates are valid unless 
the SFL splits the valence band. Otherwise much more 
charge is needed. In the first case, the density is small 
and one could expect that some surface irregularities can 
do this job. In the latter case the density is a substantial 
fraction of the defect density. We would probably notice 
some difference in the course of our experiments where 
defects appear neutral. 


E. Postulated scenario 

The above considerations strongly suggest that the 
FLP in germanium cannot be the conventional electronic 
effect adjusting the SFL to the surface state structure. 
A similar tendency has been observed for the silicon sur- 
facesiid^, but due to a larger band gap the effects are 
not that spectacular. 

It is concluded from density functional theory calcula¬ 
tions in Ref. [l^ that the Ge(OOl) surface strongly at¬ 
tracts electrons, effectively pushing the SFL towards the 
valence band. In turn, the defects differently behave in 
a near surface region compared to the bulk. This surface 
feature seems to originate from a higher amplitude of va¬ 
lence band wavefunctions close to the surface rather than 
from excessive additional surface states. In this case the 
SFL should somehow deviate upon passivation. 

The solution of this puzzle that we propose is inspired by 
putting two papers — Ref. [l^ and Ref. IIO — together. 
The first one reports ARPES experiments at higher tem¬ 
peratures on germanium and silicon samples. It is ob¬ 
served that Si(OOI) performs similar to the Ge(OOl) for 
a substantial charge accumulated in the surface states 
above the Fermi level without significant change of the 
SFL position. Note that the silicon sample under inves¬ 
tigation was flashed up to 1520 K. The latter article by 
J. L. Fitters shows that at such temperatures a depletion 
of dopant density in the subsurface layer is formed. For 
samples flashed at lower temperatures (1250 K) no mi¬ 
gration of dopants is seen. A gradient in dopant density 
is evident for the samples treated at higher temperatures. 
As a consequence, the surface Fermi level is not entirely 
determined by its bulk value and the surface state struc¬ 
ture, but also by the dopant density at the surface. The 
only discrepancy in this comparison is that the sample 
investigated in Ref. [H is phosphorous doped while in the 
latter case it is arsenic doped. 

The usual approach to the SFL pinning holds that the 
distribution of dopants is constant within the sample. 
The only way the system can respond to perturbations 
(e.g. creation of a surface) is through electronic degrees 
of freedom. As such, the system stays in a constrained. 


but not necessarily global, thermodynamic equilibrium. 
At high enough temperatures, where the diffusion is suf¬ 
ficiently quick to adjust the dopant gradient to the per¬ 
turbation, the true thermodynamic optimum is available. 
When cooling the sample the dopant gradient gets frozen. 
If the cooling is quick enough then the gradient of dopants 
reflects the equilibrium Fermi level position at the tem¬ 
perature where the diffusion stops. Such a gradient pro¬ 
file is not suited for lower temperatures making the SFL 
at lower temperatures deviate from its equilibrium posi¬ 
tion. If it does not entail charging of the surface (the case 
close to our samples) then no electrostatic field restoring 
the equilibrium SFL is created. 

Within the one-electron approximation, the Fermi level 
(chemical potential) in semiconductors is calculated by 
imposing electric neutrality conditional^ i.e. the density 
of ionized donors pdNd and holes Uh has to be equal to 
the density of ionized acceptors PaNa and electrons rig in 
conduction bands, 

PdNd + Uh = PaNa +ne, (4) 

where pd and Pa denote the relevant ionization proba¬ 
bilities and Nd and Na the dopant densities. Both the 
probabilities {pd, Pa} and densities {n^, Ue] depend on 
the Fermi level via statistical factors. As such, eq. O can 
be used to set the Fermi level position. Now, consider a 
slowly varying dopant density Nd{:x.), Na(x). Then, we 
can locally, i.e. for every x, solve the equation 

PdNdix) +nh= PaNai-X.) + Ue, (5) 

to arrive at the Fermi level as a position dependent func¬ 
tion. As such, local dopant variation impacts the electron 
Fermi level. For abrupt change of the dopant density dif¬ 
ferent calculation schemes are effective^^. 

Finally we note, that in germanium the ions are 
mobile enough to ensure swift diffusion already at 800 
K2 iiS^. This makes our hypothesis natural as we heat the 
sample to about 1000 K. Furthermore, the ion sputtering 
we use is prone to produce a considerable amount of 
defects which may be stabilized by the electrostatic field 
and change the effective doping level, see Ref. Q. 


IV. CONCLUSION 

In this paper we described electric transport experi¬ 
ments on the Ge(OOl) and Ge(001):H surfaces both in 
the classical and tunneling regimes. They show that the 
SFL is independent of the structure of surface states. 
These results cannot be understood within the conven¬ 
tional mechanism behind FLP. We argue that in our case 
a depletion of n-type dopants could be responsible for the 
experimental outcomes. 
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